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We have studied, using double ratio of QCD (spectral) sum rules, the ratio between the masses of Tec and X(3872) assuming that 
they are respectively described by the D-D* and D-D* molecular currents. We found (within our approximation) that the masses 
of these two states are almost degenerate. Since the pion exchange interaction between these mesons is exactly the same, we 
ponclude that if the observed X(3872) meson is a DD* +c.c. molecule, then the DD* molecule should also exist with approximately 
the same mass. An extension of the analysis to the fe-quark case leads to the same conclusion. We also study the SU(3) breakings 
for the Tqq/Tqq mass ratios. Motivated by the recent Belle observation of two Zt, states, we revise our determination of Xt, by 
combining results from exponential and FESR sum rules. 
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■ The existence of exotic hadrons is a long-standing problem. 
By exotic we mean a state whose quantum numbers and main 
properties cannot be explained by a simple quark-antiquark or 
three-quark configuration. 

The X(3872) resonance (assumed to be an 1^^ axial vector 
meson) has, indeed, stimulated many activities in the physics of 
hadrons. It was discovered by BELLE in B-decays UJ, and con- 
firmed by BABAR 12], CDF H and DO fl]. It is rather narrow, 
with a width < 2.3 MeV. Its most popular picture which consists 
of a molecular configuration, DD* + DD*, with J^'^ - 1+^, has 
been attributed to the narrow (< 2.3 MeV width) X(3872)[l 

' The case of the four-quark state {QQud) with quantum num- 
bers / = 0, 7 = 1 and P - +\ which, following ref.l6], we 
call Tqq, is especially interesting. As already noted previously 
||6', 7], the Tbb or Tec states with j'' - 1^ cannot split into a pair 

of two B or two D mesons which is retricted to 7^ = 0^, 2^, 

If their masses are below the BB* or DDn thresholds, these de- 
cays are also forbidden. As a result, Tqq becomes stable with 
respect to strong interaction, and must decay radiatively, or even 
weakly if the mass becomes lower than the threshold made of 
two pseudoscalar mesons. 
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'For references about some other possible interpretations of the X(3872) 
see, e.g., [5]. 



2. Tqq from potential models 

In constituent models with a flavor-independent central po- 
tential, the stability of (QQqq) configurations comes from a fa- 
vorable effect when the charge-conjugation symmetry is bro- 
ken, as noted many years ago |8]. This is the same mechanism 
by which, in QED, the loosely bound positronium molecule 
evolves into the very stable hydrogen molecule. 

It is worth noting that in the large niQ limit, the light degrees 
of freedom cannot resolve the closely bound QQ system. This 
results in bound states similar to the Aq states, with QQ playing 
the role of the heavy antiquark [9]. 

The (QQqq) states have been studied using a variety of sim- 
ple or elaborated potential models JtI [sl [lol^ 121. The corre- 
sponding four-body problem is very delicate. For instance, an 
expansion on harmonic-oscillator states was used in [iT]. It is 
eflicient for deep binding but converges very slowly for weak 
binding. If truncated, this expansion may fail to demonstrate 
stabihty with potentials that do bind, because it lacks explicit 
(Qq) - (Qq) components, which are important near threshold 
Q, and are included in the Gaussian expansion sketched in jlill 
and systematically developed in |6]. See, also, Ref. fl3| for a 
discussion about the four-quark problem. All authors agree that 
such states become bound when the quark over the antiquark 
mass ratio becomes sufliciently large. Detailed four-body cal- 
culations, using a pairwise central potential supplemented by 
a chromomagnetic interaction, indicate that Ti,h is rather well 
bound, and Tec possibly bound by a few MeV below DD*. For 
instance, the prediction of Ref. ||6[] is, in units of MeV: 

M7-„ = 3876 ~ 3905 , Mj,,,, = 10519 ~ 10651 . (1) 



A non-pairwise confinement has also been considered 01411 . in- 
spired by the large coupling regime of QCD, where it is shown 
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1 1511 that it is more favorable to build stable tetraquarks. In this 
improved quark model, as well as in conventional quark mod- 
els, it is found that (QQqq) has an energy lower than (QQqq). 

Another variant was considered in LI61I , with a chiral po- 
tential model, which includes meson-exchange forces between 
quarks, instead of the chromomagnetic interaction. 

The existence of a DD* + DD* molecule was predicted in 
ref. UtIi on the basis of the pion-exchange dynamics 0. Here, 
the pion is exchanged between the hadrons, as in the Yukawa 
theory of nuclear forces. The DD*n and DD*n vertices are 
identical, as well as the D*D*7t and D*D*n ones. There is 
only an overall change of sign, due to the G-parity of the pion. 
Therefore, if the pion-exchange dynamic^ is able to bind the 
DD* + DD* molecule, the same is true for the DD* molecule. 
The difference between these two states can only come from the 
short-range part of the interaction. 

3. Tqq from QCD (spectral) sum rules 

The first study of tetraquarks with two heavy quarks within 
QCD (spectral) sum rules (QCDSR) was done in [19] by us- 
ing diquark-antidiquark current. This study is revisited and im- 
proved in the present paper Our aim is also to compare in de- 
tail the (QQqq) and iQQq£) configurations. Such a compar- 
ison is attempted in ref. f2(J\, where the authors study heavy 
tetraquarks using a crude color-magnetic interaction, with fla- 
vor symmetry breaking corrections. They assume that the Belle 
resonance, X(3872), is a cqcq tetraquark, and use its mass as 
input to determine the mass of other tetraquark states. They 
obtain, in units of MeV: 



Mt 



3966. 



Mt 



10372 . 



(2) 



in agreement with the previous results in Eq. ([T) and the ones 
from QCD (spectral) sum rule, in units of GeV 1I19II : 



Mr,,. = 4.2 ± 0.2 , 



M- 



Thh 



10.2 ±0.3 



(3) 



The short-range part of the interaction can be tested by the QCD 
(spectral) sum rules approach 112 1142311 . Therefore, in this work, 
we study the ratio of the masses of the Tec and X(3872) states, 
by using the double ratios of sum rules (DRSR) introduced in 
1 2411 . which is widely applied for accurate determinations of 
the ratios of couplings and masses [25-31] and form factors 
1 3211 . This accuracy is reached due to partial cancellations of 
the systematics of the method and of the QCD corrections in 
the DRSR. More recently, the DRSR was used to study different 
possible currents for the X(3872) ll27ll . It was found that (within 
the accuracy of the method) the different structures (3-3 and 
6-6 tetraquarks and DD* + DD* molecule) lead to the same 
prediction for the mass. This result could indicate that the short- 
range part of the interaction alone may not be sufficient to reveal 
the nature of the X(3872). 



- For further references on this approach, see e.g. lilSll . 

^Usually, the G parity rale transforms an attractive potential into a repulsive 
one. Here, however, it only changes the sign of the transition potential DD* — > 
Z)*Z5, and thus just a phase in the two-component bound state wave function. 



3.1. Two-point functions and forms of the sum rules 

The two-point functions of the X(3872) (assumed to be an 
1^^ axial vector meson) and the Tec (assumed to be a 7^ = 1^ 
state) is defined as: 



nf(^) ^ /J d'xe''>\0\T{f.{x)jfm\0) 



= -Tiuiq'M' - ^) + n„,(?') 



(4) 



q q 

where / - X, Tec- The two invariants, EIi and llo, appearing 
in Eq. are independent and have respectively the quantum 
numbers of the spin 1 and mesons. 

We assume that the X(3872) and Tec states are described by 
the molecular currents: 

\2 



(qa(x)y5Ca(x)cb(x)y qtix)) 



eff V2 

- {qa{x)j^Ca{x)cb{x)jsqh{x)) 



and 



Jt (-^) 



A 



etf 



qaix)y5Ca(x)qh(x)yCh(x) 



(5) 



(6) 



where a and b are color indices. 

In the molecule assignement, it is assumed that there is an 
effective local current and the meson pairs are weakly bound 
by a van der Vaals force in a Fermi-like theory with a strength 
ig/^)lff which has nothing to do with the quarks and gluons 
inside each meson. 

Due to its analyticity, the correlation function. Hi, in Eq. (|4]l, 
can be written in terms of a dispersion relation: 



nu(q^) 



/-•CO 

J4»i? 



ds 



s — q^ 



(7) 



where npiis) = lm[ni,(i)] is the spectral function. 

The sum rule is obtained by evaluating the correlation func- 
tion in Eq. (HJi in two ways: using the operator product ex- 
pansion (OPE) and using the information from hadronic phe- 
nomenology. In the OPE side we work at leading order of per- 
turbation theory in a,, and we consider the contributions from 
condensates up to dimension six. In the phenomenological side, 
the correlation function is estimated by inserting intermediate 
states for the X and r„ states via their couplings A, to the molec- 
ular currents: 



<0|;;|M,) = Aie" 



(8) 



where M, = X, Tec, f are the currents in Eqs. (|5) and (|6). 

Using the ansatz: "one resonance" © "QCD continuum", 
where the QCD continuum comes from the discontinuity of the 
QCD diagrams from a continuum threshold te, the phenomeno- 
logical side of Eq. (HJl can be written as: 



nphen^ 2-. 
uv (q ) 



Mf - q^ 



(9) 



where the Lorentz structure g^y projects out the 1^ state. The 
dots denote higher axial-vector resonance contributions that 
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will be parametrized, as usual, by the QCD continuum. After 
making an inverse-Laplace (or Borel) transform on both sides, 
and transferring the continuum contribution to the QCD side, 
the moment sum rule and its ratio read: 



J 4m 



ds e Pi(s) 



4m^ 

dr 



(10) 



where r = l/M~ is the sum rule variable with M being the 
inverse-Laplace (or Borel) mass. In the following, we shall 
work with the DRSR Q: 




Mt 



My 



(11) 



3.2. QCD expression of the spectral functions 

The QCD expressions of the spectral densities of the two- 
point correlator associated to the current in Eq. (0 are given 
in ref. llssll . Up to dimension-six condensates the expressions 
associated to the current in Eq. (|6]l, in the structure g^y are: 

pis) = pP'^'is) + p^'"'\s) + 

p<'''\s)+p""\s)+p<^'>^\s), (12) 

with 



■/3)x 



Cmin J3„, 



X 



[(a + I3)m] - aPsf 



+ Ual3 + la + 5)- I5aps{\ + a + p) 



(a + /3)ml - ffySij , 
<G^>, , {g'G^} 7 da [dpi 



+513 - 3) [{a + I3)ml - aftsf + m^(l -a- pf) : 



xa^{5 + a+P) + ml(l-a-p) 
+3a + 4y6) + ap(2 + 2a+lip) + 

+ I5a-2p 



9a\2 + 



[(a+P)ml-aps]^, 



-a(l - a)s] 



/da r dp V 



12q'^ 



+ \lap- p] [{a + pWc - aps\ 



(*) = ||^ / da[nml-5a{\-a)s\, (13) 



where: nic, {g^G'), {qq), {qgcr.Gq) are respectively the charm 
quark mass, gluon condensate, light quark and mixed conden- 
sates; p indicates the violation of the four-quark vacuum satu- 
ration. The integration limits are given by: 

amin = ^(1 - V), a„ax = ^(1 + v) 

Pniin = amll{sa-ml) (14) 
where v is the c-quark velocity: 

V = -Amlls . (15) 

3.3. Tcc/X ratio of masses 

In the following, we shall extract the mass ratio Tcc/X using 
the DRSR in Eq. ( fTTT ). For the numerical analysis we shall 
introduce the renormalization group invariant quantities OTj and 
A, HQ: 

rh _ 



(-logV^A)"^^- 
{qq)(T) = -A^(-log V^Af'^^' 

(qga.GqXT) = -mlfil{-log y/^Ay'^'"' , (16) 

where Pi - -(1/2)(11 - 2n/3) is the first coefficient of the 
P function for n flavours. We have used the quark mass and 
condensate anomalous dimensions reported in 1,2 3.1 . We shall 
use the QCD parameters in Table [1] At the scale where we 
shall work, and using the parameters in Table [T] we deduce: 
p - 2.1 + 0.2, which controls the deviation from the factoriza- 
tion of the four-quark condensates. We shall not include the 
Ijq^ term discussed in 1^ 361 371. which is consistent with the LO 



approximation used here as the latter has been motivated by a 
phenomenological parametrization of the larger order terms of 
the QCD series. 

Using QCD (spectral) sum rules, one can usually estimate the 
mass of the X-meson, from the ratio "Rx in Eq- (flOl ). which is 
related to the spectral densities obtained from the current Q. A 
tetraquark current for the X(3872) was used in ref. [26]. At the 
sum rule stability point and using a slightly different (though 
consistent) set of QCD parameters than in Table [1] one obtains, 
with a good accuracy, for nic - 1 .26 GeV 1,26.] 0: 



Mv 



jUx = (3925 + 127) MeV 



(17) 



and the correlated continuum threshold value fixed simultane- 
ously by the Laplace and finite energy sum rules (FESR) sum 
rules: 



VfT^ (4.15 ±0.03) GeV. 



(18) 



"^The use of = 1.47 GeV increases the central value by about (160 ~ 200) 
MeV. 



Table 1: QCD input parameters. For the heavy quark masses, we use the range 
spanned by the running MS mass ?ng(Mg) and the on-shell mass from QCD 
(spectral) sum rules compiled in page 602,603 of the book in |23]. The values 
of A and ft^j have been obtained from a, (Mr) = 0.325(8) [ 38] and from the 
mnning masses: (m„ + mj){2) = 7.9(3) MeV |39]. The original errors have 
been multiplied by 2 for a conservative estimate of the errors. 



Parameters Values 



Ref. 



A{nf 
fig 

m, 
trie 
mi, 



'"0 



4) (324 ± 15) MeV 
(263 ± 7) MeV 
(0.114 ± 0.021) GeV 
(1.23 ~ 1.47) GeV 
(4.2 ~ 4.7) GeV 
(0.8 ± 0.2) GeV^ 
(6 ± 2) X 10-2 GeV 



138,40] 

1^3, 3£] 

[23,39,40] 

1^,39-^ 

1^,39-42] 

[28,44,45J 

|_38^4]^46-^,51.-53J 



paAddy (4.5 ± 0.3) x lO''^ GeV^ [38, 44, 46] 



In ref. Il27[l it was obtained that the DRSR for the tetraquark 
cuiTent and for the molecular current is: 



(19) 



with a negligible error. Therefore, the result in Eq. (ITSl is the 
same for the current in Eq. (|5). Although the uncertainty in 
Eq. ([TTT i is still large, considering the fact that this result was 
obtained in a Borel region where there is pole dominance and 
OPE convergence, one can say that the QCD sum rules sup- 
ports the existence of such a state and that the value obtained 
for Mx is in reasonable agreement with the experimental candi- 
date i40|]: 



Mx\exp = (3872.2 + 0.8) IVIeV 



(20) 
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Figure 1: The double ratio rj^...ix defined in Eq. 11 It as a function of r for 
= 4.15 GeV and for tvt'o values of = 1.23 (solid line) and 1.47 GeV 
(dashed line). 



We now study the DRSR of the T^dX defined in Eq. (fTTT l. 
In Fig. [1] we show the r-dependence of the ratio for yft^ = 
4.15 GeV and for two values of = 1.23 GeV and 1.47 GeV. 



^From Fig. [T]one can see that there is a r-stability around t ^ 
0.4 GeV~^ and for this value of t, we get: 

= 1.00 ±0.01 . (21) 

In Fig. |2] we show the -dependence of the ratio for t - 
0.4 GeV"^ and for two values of m,. = 1.23 GeV and 1.47 
GeV. ^From this figure one can see that the ratio increases with 
tc- However, considering the large range of tc presented in the 
figure, the ratio does not differ more than 3% from 1. 

1.05 
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0.99 



0.97 



0.95 
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Figure 2: The double ratio rr^/x as a function of for t = 0.4 GeV ^ and for 
two values of nv = 1.23 (solid line) and 1.47 GeV (dashed line). 

Our analysis has shown that the DD* + c.c. and DD* currents 
lead to the same mass predictions within the accuracy of the 
approach. The accuracy of the DRSR is bigger than the nor- 
mal QCDSR because the DRSR are less sensitive to the exact 
value and definition of the heavy quark mass and to the QCD 
continuum contributions . As mentioned before, this accuracy 
is reached due to partial cancellations of the systematics of the 
method and of the QCD corrections in the DRSR. Therefore, if 
the observed X(3872) is a molecular DD* + c.c. state its molec- 
ular cousin DD* should also be a bound state. Its mass can 
be obtained by using the experimental mass for the X(3872) in 
Eq. dlB: 

M7-„ = (3872.2 + 39.5) MeV (22) 

3.4. T hh I Xh ratio of masses 

Using the same interpolating field in Eqs. (|5]l and ^ with 
the charm quark replaced by the bottom one, we can analyse 
the DRSR: 




Mt 



(23) 



In Fig. [3] we show the r-dependence of the ratio in Eq. ( l23l l 
for -\/f7 = 10.5 GeV and for two values of m^. /,From this figure 
one can see the ratio is very stable. The same happens for the 
dependence of this ratio with tc, as can be seen by Fig.|4] We 
get: 



'•7-../x. = 1.00. 



(24) 



Therefore, we can predict the degeneracy between the masses 
of the Thb and of the Xh given in Eq. ( |25] |. 
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4. Revisiting the determination of the mass 




Figure 3: Same as Fig.[I]for n.^ix, tor = 10.5 GeV and for two values of pjgu^g 5. ;„ Q^y ^ function of r in GeV-2 from Laplace sum rule 

Mifc = 4.2 (solid line) and 4.7 GeV (dashed line). ^ 4 7 and V?;: = 1 1 GeV: long daslied (red) : (1)= perturbative 

(Pert) contribution; small dashed (blue): (2)= Pert+(c(d> + (q'5G^) contributions 
(the one of the gluon condensate is relatively negligible); continuous (olive): 
(3)=(2)+ mixed condensate (gdGd); medium dashed (black): (4)=(3)+(dd)^ . 



10.4 

tJ"(GeV) 



Figure 4: Same as Fig. |2]for rri^ijxi, forr = 0.2 GeV 
mi, = 4.2 (solid Hne) and 4.7 GeV (dashed line). 



and for two values of 
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Figure 6: in GeV as a function of in GeV from FESR (dashed curve) 
and the value at the T-stability from Laplace sum rule for mi, = 4.7 GeV. The 
OPE has been trunctaed at Z) = 5. 



The Xb was studied in ref . 02611 . At the sum rule stabiHty point 
and using the perturbative MS -mass mh{mb) - 4.24 GeV, they 
get: 

10.06 GeV < Mx„ < 10.50 GeV , (25) 

for 10.2 GeV < ^/t^ < 10.8 GeV, while combining the Laplace 
sum rule with FESR, they obtain a slightly lower but more pre- 
cise value : 



Mx„ = (10.14 ±0.10) GeV 



(26) 



We complete the previous analysis by using here the value of 
the on-shell mass m;, = 4.7 GeV due to the ambiguous defi- 
nition of the quark mass used as we work to leading order of 
radiative corrections. For a close comparison with the analysis 
in II26II . we shall work with the two-point function associated 
to the four-quark current 0. We notice that the contribution of 



^The result using the Z)S*molecule current would be the same as we have 
shown in 1 27] that the masses obtained from the four-quark and molecule cur- 
rents are degenerate. 
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the D - 6 condensate {qq}^ destabilizes the resuh [medium 
dashed (black) curve] rj which is restored if one adds the Z) = 8 
condensate given in 12611 . However, we refrain to add such a 
term due to the eventual uncertainties for controlling the high- 
dimension condensate contributions (violation of factorization, 
complete D = 8 contributions) and find safer to limit the analy- 
sis to the D - 5 contribution like in Ref. ll26ll. In this way, the 
ratio of sum rules present r-stability at about 0.1 GeV"^ (con- 
tinuous curve in Fig|5]l. We show in Fig. |6]the f^-behaviour of 
versus the continuum threshold tc, where a common solu- 
tion is obtained in units of GeV: 

Mx, = 10.50 ~ 10.78 for ^/^, = 10.5 ~ 11.0 , (27) 

which combined with the result in Eq. (l26l l leads to the conser- 
vative range of values: 



Xi, 



10.14-10.78 for Vf^ ^ My = 10.5 ~ 1 1.0 , (28) 



where one can notice the relatively small mass-difference be- 
tween Vff ™d eventually signaling the nearby location of 
the radial excitations as mentioned in 12611 . 
Very recently the Belle Collaboration studied the T(55') 
r(nS)7T^ andr(5S) ^ hb(mP)7T^ (« = 1, 2, 3 and m = 1,2) de- 
cay processes looking for resonant substructures. They found 
two narrow states in units of MeV: 



Mz, = 10610 and Mz; = 10650 , 
with a hadronic width in units of MeV: 

Tz,, = 15.6 ± 2.5 and Tz;, = 14.4 ± 3.2 
respectively 



(29) 



(30) 



channel 



I54|l. The analysis of the Zh states decay in the 
T(2S )7r^ favors the j'' - assignment, which 
is the same as the one of the X/,, although Xi, has positive 
charge conjugation and the neutral partiner of Zb should have 
negative charge conjugation. 



Considering the errors in Eq. (I28t and the small mass dif- 
ference between Xb and ^/t^, it is difficult to identify the two 
observed Zb states with the Xb. 

5. St/(3) Mass-splittings 

We extend the previous analysis to study the ratio between 
the strange T^^ {(cess)) and non-strange Tec states. 
The QCD expression for the spectral function proportional to 
nis for the current in Eq. ([6} is: 



„ OTj r c/a [ 7[ot2-q'(1 -a 



6 m^(21 -Ha-HyS) 



We have neglected the contribution of the triple gluon condensate 
8^ fabcG^) due to the (1/16^-)^ loop factor suppression compared to {qq)~. 



-6 [{a + pWc - al3s\ 



<(1 - a - /?) [(a H- P)ml - afis^ 



(31) 



We start by studying the DRSR: 




Mt 



Mt 



(32) 



We study the convergence of the DRSR versus r in Fig|2l where 
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Figure 7: Different QCD contributions to the DRSR rjsjj^^ defined in Eq. 
(32) for VfJ^ = 4.15 GeV and = 1.23 GeV: (l)=Pert+mj: dot-dashed (ma- 
roon); (2)=(1) +{qq): long-dashed (magenta); (3)=(2)+(q;jG^): dotted (blue); 
(4)=(3)+<ggG^): dashed (red); (5)=(4) +{qq}^: continuous (black). 
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Figure 8: Same as Fig.[T]for rr' /r^, defined in Eq. (32) for = 4.15 GeV 
and for two values of = 1.23 (soUd line) and 1.47 GeV (dashed line). 



the show the strength of the different contributions in the OPE. 
We note that we have r-stability around r = 0.2 GeV"^, while 
the OPE breaks down for r > 0.5 GeV 2. 
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In Fig. [8] we show the r-dependence of the ratio in Eq. ( |32] |. 
for yft^ - 4.15 GeV and for two values of nic. From this Figure 
one can deduce around the r-stability: 



= 0.95 ~ 0.98 



(33) 



which gives a smaller mass for T^^ than for Tec- This result is 
similar to the result obtained for X* in ref . (l^ . However, in this 
case, the decrease in the mass is even bigger than the obtained 
forZ^: rxs/x = 0.984 + 0.009. This resuh is consitent with what 
is obtained from the DRSR: 



Mt 



Mv 



(34) 



as can be seen by Fig. |9] 
1 
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Figure 9: Same as Fig. [T]for rjsjx'- The solid and dashed lines are for 
1.23 and 1.47 GeV respectively. 



In Fig.|9]we have used tc - 4.15 GeV. However, the result is 
very stable as a function of as can be seen in Fig.fTOl 

Since for the T^.^ state, considered as a D^D* molecule, there 
is no allowed pion exchange, one can not conclude, from the 
analysis above, that the 7"*^. should be more deeply bound than 
the Tec- On the contrary, if the pion exchange is important for 
binding the two mesons, the may not be bound. 

6. Conclusion 

In conclusion, we have studied the mass of the Tec using dou- 
ble ratios of sum rules (DRSR), which are more accurate than 
the usual simple ratios used in the literature. We found that 
the molecular currents DD* + c.c. and DD* lead to (almost) the 
same mass predictions within the accuracy of the method. Since 
the pion exchange interaction between these mesons is exactly 
the same, we conclude that if the observed X(3872) meson is a 
DD* + c.c. molecule, then the DD* molecule should also exist 
with approximately the same mass. A recent estimate of the 
production rate indicates that these states could be seen in LHC 
experiments 1155 1. 
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0.92 



4.4 4.6 

t„"'(GeV) 

Figure 10: Same as Fig. |2]for rjs jx^- The solid and dashed lines are for : 
1.23 and 1.47 GeV respectively. 



We have also studied the double ratio rj^^ixt using molecular 
currents BB* and BB* + c.c. for Ti,b and Xt respectively. In 
this case the degeneracy between the two masses is even better 
than in the charm case. Therefore, we also conclude for the 
bottom case that if a molecular state BB* + c.c. exist, then the 
BB* molecule should also exist with the same mass. 
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